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The leading term of the asymptotic of quasinormal modes in the Schwarzschild 
background, iv n = — in/2, is obtained in two straightforward analytical ways for 
■ arbitrary spins. One of these approaches requires almost no calculations. As simply 

we demonstrate that for any odd integer spin, described by the Teukolsky equa- 
^sO | tion, the first correction to the leading term vanishes. Then, this correction for 

half- integer spins is obtained in a slightly more intricate way. At last, we derive 
analytically the general expression for the first correction for all spins, described by 
^sD ' the Teukolsky equation. 

in : 
o ■ 

1 Introduction 

in : 

The investigation of perturbations of various fields in the Schwarzschild background was 
started in |UI2]- Quasinormal modes (QNM) are the eigenmodes of the homogeneous wave 
equations, describing these perturbations, with the boundary conditions corresponding to 
outgoing waves at the spatial infinity and incoming waves at the horizon. The interest to 
QNMs was initiated by HE]. 

Two boundary conditions make the frequency spectrum uj n of QNMs discrete. The 
asymptotic form of this spectrum for gravitational and scalar perturbations of the Schwarz- 



schild background was found at first numerically in jSJ E] : 

u n = -- (n+ -) + 0.087424, n -»• oo , s = 0, 2 . (1) 



2 V 2, 

Here and below the gravitational radius r g is put to unity; s is the spin of the perturbation. 
This result up to now serves as a touch stone for investigations in the field. 

A curious observation was made in |7J: the real constant in (JTJ can be presented as 

Reu n = ^ = T H In 3, (2) 

where Tjj is the Hawking temperature (T# = l/(87ckM) in the common units) 3 . Then, 
expression (j2J for the asymptotic of Reu n was derived in [0] by solving approximately the 
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3 It was also conjectured in [J] that the asymptotic value J5J) for Rew„ is of a crucial importance for 
the quantization of gravitational field, fixing the value of the so-called Barbero - Immirzi parameter. In 
spite of being very popular, this idea is not in fact dictated by any sound physical arguments; quite the 
contrary, it is in conflict with them [Sj. 
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recursion relations used previously in the numerical calculations. In the next paper JUj 
formula (J2J) was derived analytically. Besides, in jS] the following result was obtained for 
the asymptotic of QNMs for spin 1: 

i 

uj n = —-n, Rea; n ^0, n — > oo , s — 1. (3) 

Asymptotic (JHJ) was also obtained numerically in [IT] . 

While the results (pQ) and (JHJ) for integer s are firmly established now, it is not the 
case for spin 1/2. Two different approaches 4 used in |T2] result in the interval two times 
smaller than those for integer spins, namely: 

uj n = — -n, n — > oo , s = 1/2. (4) 

On the other hand, numerical calculations in |HJ result in spectrum 

u n = — -n, Reu n — >0, n — > oo , s = 1/2. (5) 

One of the motivations of our work was the resolution of this discrepancy; we not only 
confirm below equation but find also first nonvanishing correction to it. 

We consider the QNM problem in various analytical approaches. Two of them, rather 
simple and straightforward, give in fact only the leading asymptotic, u n = — in/2 for any 
spin. In the third approach, based on the Teukolsky equation, we at first demonstrate 
as easily that equation (J3J) is accurate for arbitrary odd integer spins, i.e. that first 
subleading correction to it vanishes. Then, with somewhat more efforts, we obtain this 
correction for half-integer spins. At last, we derive, in a more involved way, the unified 
general expression for the next term in the asymptotic values of the QNMs for all spins, 
which was conjectured previously in [TU] . 



2 Quasinormal modes in Regge — Wheeler formalism 

The Regge - Wheeler equation for the radial function \l/ corresponding to the angular 
momentum j of a field with integer spin s (s = 0, 1, 2 ; j > s) is written usually as 



d 2 ^ 



1 



* = 



(6) 



Its analogue for s — 1/2 (again the angular momentum j > s), written for the standard 
representation of the Dirac 7-matrices and states of definite parity, is 




3/2' 



r J 



* = 0; (7) 



4 We believe that one of them, despite being rather popular, can be dismissed at once. It is based on 
the analysis of the location of the poles of the scattering amplitude, which by itself causes no objections. 
However, following Ji] _ CHjj the authors of ^2] analyze the poles of the corresponding Born amplitude. 
Meanwhile, the Born approximation by itself implies that the amplitude of the scattered wave is small. 
Therefore, its poles have no real meaning. Any coincidence between their position and that of the poles 
of a true amplitude is an accident only. 
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here k = ±(j + 1/2), with the sign depending on the parity of the state considered (this 
sign is irrelevant for our problem). The presence of the terms with fractional powers of r 
and r — 1 in equation (|7J) is quite natural since wave equations for half-integer spins are 
written via tetrads which are roughly square roots of metric. 5 

In both equations, (JUJ) and J7J), r is treated as a function of the so-called "tortoise" 
coordinate z. They are related as follows: z = r + ln(r — 1), so that z — > oo for r — > oo, 
and z — > — oo for r — > 1. The boundary conditions for QNMs of © and (J7J) are 



$?(z) ~ e 



±00. 



Here, for our purpose, it is convenient to go over in both equations, © and (J7|), to 
the usual coordinate r and to the new radial function u(r) related to \1/ as follows: 



,1/2 



(9) 



The obtained equations for u{r) can be rewritten as 




(10) 



+ 



3 + 



3 1 

4 j*2 



1 



ft 







1/2 



r 3/2^ r _ x) 1 / 2 2 r 3 / 2 (r — l) 3 / 2 

We are interested in the solutions of equations (|r0*j) and (fTTjl in the interval 1 < r < oo 
for | a; | — > oo. Obviously, all the terms singular at r — > 0, in both these equations, are 
relatively small in this interval if | u\ — > oo. 6 Therefore, these terms can be safely omitted, 
and we arrive at the following universal truncated wave equation for all spins: 



d 2 u 
dr' 1 



+ 



co 2 + 



2u 2 oo 2 + 1/4 



+ 



(r 



u 



0. 



(12) 



We have omitted here also the terms —(J + 1/2) 2 + s 2 and —(J + 1/2) 2 + 1/2 in the 
coefficients at l/(r — 1) in (fTUj) and (fTTj) . respectively. Though these terms could be easily 

5 We mention here another rather popular, but false belief, namely, that equation © applies to half- 
integer s as well. The explicit difference between (0 and (0) demonstrates that this idea is wrong. 



6 In particular, in equation (flip 



and | k | r- 3/2 {r - iy 3/2 < | u 2 \ (r - 1)" 



for the interval 1 < r < oo . 
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included into the solutions, they would result in corrections to Im uo n on the order of 1/n 
only, which are negligible as compared to the leading term ~ n. 

We retain however the term 1/4 in the coefficient at l/(r — l) 2 in (|12j). Otherwise 
the wave function asymptotic for z — > — oo would be e ~ luJZ + 1 / 2 ^ instead of e~ tu)Z . In other 
words, the effective potential in the initial Regge - Wheeler equations would not 
vanish for z — > — oo, but would tend instead to 1/4. Indeed, the wave function asymptotic 
for z — > — oo is determined by the discussed coefficient at l/(r— l) 2 . Since the coefficient 
oj 2 + 1/4 at l/(r — l) 2 in (|T2j) corresponds to J 2 in equations ©, ©, then obviously the 
coefficient u 2 in (fT2|) would correspond to u 2 — 1/4 in ©, (|IJ). 

To summarize, it is only natural that equation ()12j) . essentially semiclassical one (due 
to the assumption \ uo\ ^> 1), is universal, i.e. independent of spin s. Moreover, even if 
one assumes that j ^> 1 as well (i.e. gives up the condition j <^ \ u\ used in (|T2"j) ). the 
resulting, again semiclassical equation 

d 2 u f 2 2lu 2 1 / 2 1\ If 1\ 2 } . , 

+ ^ + (^r + 4hK^l j+ 2j j (13) 

is still universal, i.e. spin-independent. 

We address now the eigenvalues of equation (fl2"j). Its two independent solutions can 
be conveniently expressed via the Whittaker functions W\ tfl (x) (see, e.g., [Hj). They are 

W^{-2 1 cu(r - 1)) , W^U2 1 u{r - 1)) . 

With their different asymptotic for r — > oo, 

W^-2iuj{r - 1)) - e ^r + Hr-i)] = 

W^{2iu{r - 1)) -> e -*-KMr-i)l = e — ; 

these solutions are obviously independent. On the other hand, the second one does not 
comply with boundary condition (|SJ) and therefore should be excluded. 
As to the first solution, its limit for r — > 1 is 

W iuij U-2iu{r-l)) — > 



r(l/2-2iw) T(l/2) 1 v ;J v ; 

When going over to the function \l/ used in the "tortoise" coordinate z (see (jHJ)), the 
overall factor (r — 1) 1//2 in this expression cancels, and (r — l) ±4a; goes over into e ±luz for 
r — > 1. To comply with the boundary condition at the horizon, one should get rid of the 
first term in equation (JHJ). To this end, recalling that T(—n) has poles for integer positive 
n, we put 1/2 — 2iu = —n, or u n = — (i/2) (n + 1/2). 

In fact, equation (|12|) by itself was obtained from (jlUJl and (|llj) under the assumption 
|c<j n | — > oo, or n > 1. Therefore, in this way we can guarantee, for the initial problem, 
only that 

i 

u> n = --n, n> 1, (15) 

for all spins. 
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Figure 1: Singular point of equation (JT2J), cut, and closed contour 

Though less accurate than quantization rules (JIJ and Q, this one is still quite sufficient 
for insisting that the correct leading term in the quantization rule for spin 1/2 is (0), but 
not ©. 

In conclusion of this section, we demonstrate, with a relatively simple example of 
truncated equation an analytical approach, that will allow later, for more accurate 
treatment of the wave equations, to find not only the leading term ~ n in the asymptotic 
of QNMs, but as well the next, constant one. The method goes back to [TUj where it was 
applied to the Regge - Wheeler equation for s — 0, 2. After finding in the present section 
by this method the eigenvalues of equation (fT2j) . we will apply below the technique to the 
Teukolsky equation for arbitrary spins. Our line of reasoning differs from that of [TP] . 

Equation (fT2*|) has two singular points, r = 1 and r = oo. We connect them by a cut in 
the complex plane r going, for instance, from r = 1 along the real axis to the right (solid 
line in Fig. 1). Let us consider the closed contour marked by the dashed line in Fig. 1. 
Since there is no singularity inside it, the solution at some point on this contour, after 
going around the contour, comes back to its initial value, which means that the phase of 
this solution changes by 2 7r n, n = 0, ±1, ±2, 

When we follow an arc of a large radius r ^> 1, where the asymptotic solution is 
e lujr r luJ ^ i. e . go around the singular point at infinity, the wave function acquires the phase 
S(oo) = 2niui. 

Then we go around the branch point r = 1 by following an arc of a small radius. 
Here, due to the asymptotic solution v (r) = (r — l)-^ +1 / 2 , the wave function acquires 
the phase 6(1) = 2n(iu — 1/2). As to the paths along the cut, they generate no phase 
at all. Indeed, since r = 1 is a regular singular point, the wave function can be written 
as u[r) = v(r)w(r), where w(r) is analytic at r = 1. The phase of v (r) = (r — i)-' luJ + 1 / 2 
remains constant along the paths adjacent to the cut, as well as the phase of r — 1. As to 
the analytic function w(r), it obviously cannot acquire any phase after going around the 



cut. In other words, effectively for our purpose, the regular singular point r = 1 behaves 
as if it were an isolated singularity. 

Thus, going counter-clockwise around the considered closed contour in the complex 
plane, one obtains 

8(00) + 8(1) = Airiuj — 7T = 2im , 

or the quantization rule 

0J n = - ~ (n + 1/2) . 

Being interested in the solutions decreasing in time, we choose here positive n (and of 
course large ones). Again, one can guarantee here the leading term only, u n = — in/2. 



3 Teukolsky equation. Quasinormal modes 
of odd integer spin 

Now we address the problem of the next, subleading correction, of zeroth order in n, to 
formula (fT5j). It is only natural to expect that this correction is spin-dependent. So, to 
investigate it we will use the Teukolsky equation. As distinct from the Regge - Wheeler 
equation, this one describes in a unified way both integer and half-integer spins, ranging 
at least from s = to s = 2 [H]-|2n|- Previously, the Teukolsky equation was used in jS] 
for numerical calculations of QNMs. 

In the Schwarzschild background the Teukolsky equation for a massless field is 

A^ + (l- a )(2r-l)^ + tf(r)J2 = 0, (16) 

where 

a / \ , „ N TT , x - r(2r - 3)iuj s + r 3 uj 2 ... 
A(r)=r(r-1), U(r) = ^ A js , A js = (j + s)(j - s + 1) . 

Obviously, for a given spin s the QNMs are independent of helicity. 

With the tortoise coordinate z(r) = r + ln(r — 1) and new function x( r ) = rA~ s / 2 R(r), 
one obtains the following standard form for this equation: 



fx 

dz 2 



+ [u 2 -V(r)] X = 0, (17) 



with the effective potential 



s 2 -A A 2 s -s + s 2 -l A 2 js -s + s 2 -3iu;s 2iuos 
V ^=^i ' -3 + ~ ~2 +^T- • ( 18 ) 

Clearly, for s = V(r) is real, and equation ()17j) coincides with the scalar version 
of the Regge - Wheeler equation (jOJ). On the other hand, the Teukolsky equation for 
s = 1/2 coincides with the second order equation for a massless Dirac field in the chiral 
representation; of course, the latter differs from equation ((7J) written in the standard 
representation. 
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Figure 2: Singular points of equation (|T7jl . cuts, and closed contour 

The asymptotic behavior of QNMs of the Teukolsky equation is 

x ~A s/2 e ±iuJZ for ujz -> ±00. (19) 

In principle, here the idea of calculating the eigenmodes will be the same as the second 
one above, in the case of truncated equation (fT2"|) . We choose a closed contour without 
any singularity inside it (dashed line in Fig. 2)), calculate the phase of the wave function 
acquired after going around the contour, and equate this phase to 27m, n = 0, ±1, ±2, ... . 
However, this problem for the Teukolsky equation is in general much more involved than 
that for the truncated equation (|T2*|) since equation (JTTj) has three singular points: r = 0, 
r = 1, and r = 00. As previously, we choose a cut in the complex plane r going from r = 1 
to r = 00 along the real axis to the right (solid line in Fig. 2). As to the cut starting at 
r = (another solid line in Fig. 2), the choice of its location will be discussed later. 7 

The treatment of the first cut is practically the same as in the simple case of equation 
(112}) . The present boundary condition (|19|) means that the asymptotic solution for r — > 1 
looks here as (r — l)- ?w + s / 2 . Correspondingly, when going around the branch point r = 1 
along an arc of a small radius, the wave function acquires the phase 

5(1) = 2<k(iuj - s/2) . (20) 

In this well, the regular singular point r = 1 is effectively equivalent, for our 

purpose, to an isolated singularity. 

Situation with the cut starting at r = is more complicated. The problem is that we 
have no a priory boundary condition at r = 0, but still have to find the phase acquired 
when going around this cut. 

However, for s — 1 (and for any other odd integer spin that is described by the Teukol- 
sky equation) the discussed correction can be found easily. For r<l, two independent 

7 Dotted lines in Fig. 2 are the level lines, see below. 
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solutions of equation (fTTj) with potential (fTHj) are 



X+ 



+ + s ' 2 and 



X- 



r 



l-s/2 



(21) 



Since r = is a regular singular point of the Teukolsky equation, the exact general solution 
of this equation can be presented as follows: 



x(r)=r W2 



a * r ' k + rS J2 b 



_k=0 



k=0 



-s/2 



kT 



fe+s+1 



k=0 



k=0 



(22) 



With an odd integer s, the singularity of this solution at r = is due to the overall factor 
r~ s / 2 , or to r s l 2 if by some reasons X- vanishes (we will see in the next section that just 
this is the case) . Correspondingly, the phase acquired by solution (|2Tj) as a result of going 
around the branch point r = is 



5(0) = ±7TS . 



(23) 



In fact, the sign in this expression does not matter for our problem since the two options 
differ by 2ns with an integer s. And if necessary, one can always shift the initial n 3> 1 
in the quantization rule by an integer s. 

At last, the asymptotic solution at infinity 



X~r s+ "e ,ur , r^oo, (24) 

after going by 2ir around the arc of infinitely large radius, acquires the phase 

5{oo) = 2tt(zo; + s) . (25) 

With the phase (|2Uj) generated by the branch point r = 1, the total result of going around 
the closed contour is 

5(1) + 5(oo) + 5(0) =4niu 



(we have chosen here the sign minus in (|23))). Equating this expression, as above, with 
27m, we arrive at the quantization rule already mentioned in Introduction (see 0), but 
now for any odd spin: 

% 

oj n = n, s = l,3, .... (26) 



4 Quasinormal modes of half-integer spin 

The situation is somewhat more complicated for half-integer spins. Here we have to find 
out which of the two solutions that behave at r — ► as r 1+s / 2 and r 1-s / 2 , respectively, is 
the true one. To this end we need to match the solution for | r| <C 1 to that for \r\ 3> 1. 
Fortunately, in the limit | u\ ^> 1 it can be done analytically (here we follow the idea used 
in JU] for finding the QNMs of Regge - Wheeler equation (|B|)). 

Here and below, to investigate the singularity at r = 0, it is convenient to shift 
z —>■ z + iir, so that now z(r) = r + ln(l — r), and in the limit r C 1 we have z(r) = —r 2 /2. 
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Introducing new variable p = uz{r) = — ur 2 /2, we transform equation (j!8|) in the limit 

I U)\ > 1 to 8 



dx 



2 



dp 



+ 



3is 4 — s 



+ 



2p 16 p 



2 



X = 0. (27) 



Independent solutions of equation (J2*7)l are the Whittaker functions 

W^A-2ip) and W_as A2ip) . (28) 

Though derived for | r\ <C 1, these solutions are valid also for \p\ = | ur 2 /2 | 3> 1, if \w\ 
is sufficiently large. Their asymptotic behavior for \p\ = \uz\ — > oo is, respectively, 

p^e ip and p-* e -v\ (29) 

To choose the appropriate solution we compare the asymptotic behavior (|29|) with that 
of the solution of exact equation (|17|) for \ujz\ — > oo, as given in (|19|). With the leading 
asymptotic for QNMs already established in section 2 for arbitrary spins, uo n ~ — in/2, 
the solution (|T9"j) is exponentially small in the left half-plane r. Therefore, we have to 
choose here the exponentially small solution of equation (|2*7jl . In this way we arrive at 



X{p) = Wto t .(-2ip) ~ p^e ip . (30) 



We note that matching of the two solutions, (fT9~j) and (J5U|) . is not precluded by the 
fact that the pre-exponential factors in their asymptotic, ~ [ujz) s and p^ ~ (qjz)~£, 
respectively, are different. This difference is only natural since the factor (ujz) s is due to the 
term 2ius/r in equation (|18j) . and the factor {ojz)^t is due to the term —3is/2p = 3is/ur 2 
in equation (|27|) . The coincidence of the exponentials themselves is quite sufficient reason 
to believe that it is just ()3U|) that reproduces the behavior of the exact solution for |r| — > 0. 

The functions W^^y) with a given asymptotic behavior for \y\ — > oo can be expressed 
via other linearly independent solutions M\ )fl of the Whittaker equation with a definite 
behavior for \y\ — > 0. These solutions are [Tj\ 

M x> „(y) = y^ e-tf $(1/2 + - A, 1 + 2p, y) , (31) 

where $(a,b,y) is the confluent hypergeometric function. Functions W\ tli and Ma )M are 
related as follows [T7j : 

^,(y) = r(1/ r 2 ( :^_ A) M,,(y) + T{1 ™_ A) M,-.fa) . (32) 

In the present case we have 

W^(-2ip) = ^0^M^(-^p) + r(1 ^ /2) M ¥ ,_.(-2ip) . (33) 



8 If one retains one more term in the expansion of z(r), i.e. with z{r) = — r 2 /2 — r 2 /3 , a correction 

l-a 2 + 3j(j + l) 

ev^p 3 / 2 ^ 1 / 2 ' 

of first order in |lj| _1 / 2 , arises in square brackets of equation (|27J) . It coincides with the corresponding 
perturbation obtained in [21] that generates corrections ~ n -1 / 2 to u> n for s = 0, 2 in the Regge - Wheeler 
formalism. 
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For half- integer s, the first term in this expression vanishes since V(z) turns to infinity 
for negative integer z. So, here our solution (J3UJ), with the account for ([3~T]) and with 
p ~ r 2 , behaves for r — * as 

X ~ r X - s/2 • (34) 
Here again (see (j23|) ) the phase due to the branch point r = is 

5(0) = tcs, (35) 

but now s is half-integer. We have again the same 5(1) and 5(oo) as those for odd spin 
(see (I2UJ) and (|25)1. respectively), and the same quantization condition 

5(1) + 5(oo) + 5(0) = 27m. 

At last, shifting the initial n by the integer part of (now half-integer) s, we arrive at the 
quantization rule for any half-integer spin: 

w »=-^( n+ 0> ^ = 1/2,3/2,.... (36) 

Of course, equation (fHHJ) can be directly employed also for odd integer s. In this case 
the second term in this equation vanishes, and now obvious line of reasoning results in 
formula (J2HJ). 



5 Quasinormal modes of arbitrary spin 

The problem for arbitrary s, including the case of direct physical interest, that of s = 
0, 2, requires more sophisticated approach. In this general case both terms in the rhs of 
equation (J3*3*j) survive, so that the solution near the origin contains both powers of r: 

x = ar 1 + s/2 + 6r 1 - s/2 . (37) 

Obviously, the rotation of this expression by 2ir around the branch point r = in no way 
can result in its multiplication by some factor, i.e. the limit (}3*Tj) of the solution for small 
r cannot transform into itself under this procedure. 9 However, such a transformation 
does exist for the solution far away from the origin. The rotation gets possible due to the 
Stokes phenomenon, rather well-known in mathematical physics. 

We will not discuss this phenomenon in general, but instead will demonstrate directly 
how it works, by solving our problem. We consider the approximate solution Wss s(—2ip) 

4 ' 4 

for generic s and p (though confine of course to small r). The analytic continuation to 
the specific values of spin will be performed only in the final result which is a smooth 
function of s. 

At first we discuss the position of those lines in the complex plane r where Im p = 
Im (ojz) = 0. These four level lines in the complex r plane are presented in Fig. 2 

9 An additional problem arises for even s. In this case T(— s/2) in equation (|33|l turns to infinity (and 
for s = 0, r(s/2) turns to infinity as well), so that this solution stays finite due to a delicate cancellation 
between two terms in the rhs of equation Therefore, to obtain the explicit form of the solution for 

even s one should perform a careful limiting transition. Anyway, this solution for small r also does not 
transform into itself under the rotation around the origin. 
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(dotted lines therein). Their behavior corresponds to the leading asymptotic for QNMs, 
already established in section 2 for arbitrary spins: uo n ~ —in/2. For small | r| (still, 
\p\ = | ur 2 /2 | can be large!) the definition of level lines in the complex r plane reduces 
to Im(u;r 2 ) = 0, or argr = — 1/2 aigu. Two of them, c and d, going to the right from the 
origin, with argr = ±7r/4, are of no special interest to us. We will be interested mainly 
in the two level lines, a and b, that become vertical at infinity. In the sector between 
these lines the exact solution is exponentially small for r — > oo (see (I24|) ) . 

Our solution ()30|) has a cut going from the branch point r = to infinity. This cut 
should be chosen in a judicious, self-consistent way. 

For instance, it cannot be drawn in the sector where the solution is exponentially 
small. Indeed, as it was demonstrated above, if we went around such a cut starting from 
the small solution (JHUj) . we would arrive in the result at a linear combination of both small 
and large solutions. But a large solution should not exist in this sector. 

By the same reason, if starting from the real positive r axis we go in the positive 
direction, i.e. counter-clockwise, along the contour of large r, the cut cannot be drawn 
along the level line a or in the sector to the right of it. 

Neither, with this direction, should we draw the cut in the sector between the level 
lines b and c where the solution is exponentially large at infinity. In this sector we cannot 
guarantee the absence of an exponentially small admixture to the right of the level line b; 
then, after going around the cut, the correct solution would be completely distorted. 

Thus, for the counter-clockwise direction, the only consistent choice for the cut is that 
along the level line b. Just in this way we will proceed. 10 

So, let us go from the real r axis in the counter-clockwise direction along a contour 
of large r. We reach the level line b, and then proceed along its upper side. At a small 
distance from the origin, we follow an arc of a radius r <C 1. Then we come back along 
the lower side of the cut to the arc of a large radius r ^> 1. At last, we close the contour 
by going along this arc, and then around the cut starting at r = 1. 

As to the cut starting from the origin, here again the contributions to the acquired 
phase from the upper and lower sides of the cut cancel. So, we have to find only the 
phase generated by the rotation around the branch point r = 0. The rotation angle 
here is — 2tt in the r plane, which corresponds to — Att in the p plane. To calculate the 
mentioned phase, we use solution (|3U|) in the limit p ^> 1, i.e. we work in the interval 
| uj |~ 1/2 < r < 1 (recall that | uj | > 1). 

To perform the rotation by 27r in the r-plane, or by 4tt in the p-plane, we note first of 
all that functions transform under the rotations in a simple way. Indeed, according 
to equation (j3T|) . 

M x ^{ye-^)=e- 4 ^M K ,{y). 

Then we need the relation inverse to (|3~2j). to express back Ma i(U via W\^. Its form 
depends on aigy (see [T7j, 9.233.1, 9.233.2). In our case, the required initial value of 
y = —2ip corresponds to the upper side of the cut along the level line b. Since arg?/ 
remains constant along a level line, it can be found most easily for |r| — > oo where on this 
line argr = 3n/2. In such a way, we have here 

71 71 

argy = arg (—2ip) = — — + argcu + arg r = — . 

10 Quite analogous arguments demonstrate that if one goes from the real r axis in the negative, clockwise 
direction, the cut should be chosen along the level line o. 
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With this argy, the inverse relation reads 

^ r(1 7 2 1 + + r + ) A ) -"" /2+ ^' y ^'^ -£<«*«<£■ ^ 

Now we go back in the rotated solution 

+ r(i/2 2 + ^x) ei '*'' M ^- 2 ^ (39) 

to functions W±\ t ±^(±2ip) by means of ()38|) . To simplify the result of this transformation, 
we note first of all that the contributions to the result originating from the first term in 
equation (jHEJ) (containing W-\ tfl (e~ m y)), are exponentially small along all the path leading 
from the line b to the real axis; besides they vanish of course on this axis. So, these terms 
can be neglected at all. 11 Then, we need the result only in the limit \p\ ^> 1, where 
W\ )fi (y) = W\ t -[j,(y). In this way, we arrive at relation 

x (pe- iiT ) = - e- ins (1 + 2 cos its) W M i(-2tp) . (40) 

4 ' 4 

The coefficient — e~ l7TS (l + 2 cos7rs) here results from trivial, but rather tedious transfor- 
mations with T-functions, sines, and cosines. This coefficient can be rewritten as e % 5 ^ , 
where 

S(0) = 7T - 7TS - iln(l + 2C0S7TS) (41) 

is the phase acquired by the solution when following the arc of a small radius r < 1 
around the origin r = 0. 

As usual, the quantization condition for u n is 

6(1) + 6(oo) + 5(0) = 2vm. (42) 

Finally, with ()20|). (|23j) . and (|41j). we obtain analytically the universal formula 

1 ( 1\ 1 , . 

v n =—n[n+-)+ — m(l + 2 cos 71 s) , n —> oo (43) 

2 \ 2 1 An 



for eigenmodes of any spin s described by the Teukolsky equation. 
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1:L In the region where \y\ ^> 1 and the asymptotic form of the function is used, we in fact have neglected 
already small power-like corrections to this form. So much the more, we can and even should neglect 
exponentially small corrections to it. This is the Stokes phenomenon at work. 

12 Of course, the same result arises when going in the opposite direction along the contour in the 
complex r plane, with a cut made along the level line a. However, in this case equation (I38|) modifies to: 

M, ( v ) - F ( 1 + 2/i ) ur x (e i7T v) + T ( 1 + 2 ^ e -«r(l/2+M-A) w / \ . 

Mx,^y) r(1/2 + /i _ A ) e W-A^te y)+ r(1/2 + /i + A) e W\,AV), 

37T 7T 

— 2~ < ar §2/ < 2 ■ 



12 



6 Conclusions 



Now few comments on our results. 

It is not clear whether the Teukolsky equation is valid for s > 2 . 

For s = 0, 2 formula (|4Hj) was derived previously in (HI E]j in the Regge - Wheeler 
formalism. For these spins it gives 

if l \ 1 i 
w n =--ln+-l+ — ln3, n-^oo. 

The result 

z 

uj n = - -n , n — > oo , 

for s = 1 was previously obtained in [HIE]- It is derived in an elementary way in section 3, 
confirmed in section 4, and follows immediately from ([43)1 . 

For half- integer spins, simple calculation in section 4, as well as formula (|43|). give 

cg>„ = — — ( n + — ) , n — > oo , s = 1/2,3/2. 



2 V 2, 

For s = 1/2 it not only confirms the conclusion of ^3], thus resolving the controversy on 
the matter, but contains also first nonvanishing correction to the leading term; the result 
for this correction is new. Our total result for s = 3/2 is also new. Quite recently, both 
results were confirmed in 
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